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Abstract 

It  is  shown  that  a  solution  in  closed  form  can  be  obtained  for  the  problem 
of  finding  the  pressure  distribution  on  a  circular  disc  in  an  arbitrary  axi- 
symmetrical  sound  wave  when  it  is  assumed  that  the  scattered  field  behaves  like 
a  standing  wave  at  infinity  and  approaches  the  static  limit  as  the  wavelength 
becomes  infinitely  long.  This  result  is  used  to  convert  the  standard  integral 
equation  of  the  first  kind,  which  arises  when  the  scattered  field  behaves  like  a 
radiating  wave,  into  an  integral  equation  of  the  second  kind. 

The  new  integral  equation  has  a  particularly  simple  kernel  so  that,  although 
it  is  difficult  to  solve  exactly,  it  can  be  solved  approximately  with  ease.  It  is 
found  that  one  approximation  leads  to  an  expression  (in  terms  of  well-tabulated 
functions)  for  the  scattering  coefficient  which  agrees  adequately  with  the  values 
obtained  from  the  spheroidal  function  expansion  provided  that  the  product  of  the 
wave  number  and  radius  of  the  disc  is  less  than  3.  The  approximation  is  capable 
of  systematic  extension  to  cover  larger  values  of  the  product  of  wave  number  and 
radius.  The  approximation  is  also  used  to  determine  the  scattering  coefficient 
when  the  incident  field  is  derived  from  a  point  source  on  the  axis  of  the  disc. 

A  general  theory  is  developed  of  the  corresponding  integral  equation  which 
would  arise  in  aiiy  scattering  problem.  It  is  shown  that  upper  bounds  can  be  ob- 
tained for  the  scattered  amplitude  and  the  scattering  coefficient.  Numerical 
values  are  given  for  a  circular  disc  in  a  normally  incident  plane  sound  wave  and 
in  the  field  due  to  an  axial  point  source. 

Further,  it  is  shown  that  it  is  always  possible  to  estimate  the  error  made 
in  finding  the  source  distribution  by  solving  the  integral  equation  approximately. 
From  this,  estimates  of  the  error  in  the  scattered  amplitude  and  ecattering  coeffi- 
cient are  deduced.  Numerical  results  are  given  for  a  circular  disc  in  a  normally 
incident  plane  sound  wave. 
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Introduction 

The  problem  of  the  diffraction  by  a  circular  disc  or  aperture  has  attracted 
much  attention  in  recent  years.     It  is  the  simplest  problem  in  which  there  if  a 
diffracting  edge  of  finite  length  and  therefore  serves  as  a  convenient  test 
for  trying  out  methods  which  may  be   applicable  to  more  complicated  problems. 
These  methods  have  been  reviewed  by  Bouwkampl-  -»,   so  that  here  it  is  necessary 
only  to  describe  briefly  one  or  two  of  them.     The  direct  method  of  separation  of 
variables  leads  a  solution  in  spheroidal  functions  ■iriiich  was  given  by  Bouwkamp  ^  -I , 
In  principle,  all  that  is  required  can  be  obtained  provided  that  sufficient  com- 
putation and  information  about  spheroidal  functions  is  available.     To  date,  numer- 
ical results  have  been  given  in  the  range  0  <  a  <  10,  where  a  is  the  product  of  the 
wave  number  and  the   radius  of  the  disc.     Thus,  in  practice,  the  method  of  separation 
of  variables  is  most  useful  at  the  low-frequency  end  of  the   frequency  scale.     It  is, 
however,   applicable  only  to   special  shapes.     In  an  attempt  to  provide  a  method  which 
was  not  restricted  to  shape,  Le vine  and  Schwinger  ^-=^-1  introduced  a  variational  princi- 
ple.    They  were  able  to  obtain  good  agreement  with  Bouwkarap,  but  it  has  been  shown 
(Jones'--')  that,  in  general,  only  a  special  choice  of  trial  function  (which  coin- 
cides with  Levine  and  Schwinger's  choice  in  the  problem  of  the   circle)  can  lead 
to  accxjracy  (in  a  certain  sense)   at  low  frequencies.     The  determination  of  the 
special  trial  function  is  equivalent  to  solving  the  problem  by  means  of  a  power 
series  expansion  in  terms  of  the  parameter  a.     The  expansion  correct  to  terms  of 
0(a  )  was  derived  by  BouwkampL  K     The  range   of  validity  of  such  an  expansion  is 
not  known,   although  there  is  some  indication  from  the  following  work  and  that  of 
Magnus  *■  -■  that  it  does  not  extend  much  beyond  a  "  1,     Another  method,  using  an 
assumed  aperture  distribution  in  terms  of  hypergeometric  polynomials,  has  been 
discussed  by  Chako  L -I   and  by  Nomura  and  Katsura'-  -1.     The  application  of  the  Wiener- 
Hopf  technique  has  been  considered  by  Leitner  and  Wells  •-  J  and  leads  to  mathemati- 
cal difficulties  similar  to  those  encountered  in  diffraction  by  a  finite  waveguide 

[91 
Jones'-  -'• 
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At  high  frequencies  there  is  no  exact  mathematical  sclutlcn  irtiich  gives  results 
correct  to  O.l/j  (say)  in  practice,  although  in  principle  these  could  be  obtained 
froB  the  spheroidal  function  solution.     There  are   available,  however,  approximations 
based  on  physical  reasoning,     "nie  classical  method  is  due  to  Kirchhoff,  but  more 
convincing  explanations  of  what  happens  have  been  given  by  ClemBOw   •-    ^   and  Keller  L"^-", 
Their  explanations  have  the  virtues  of  being  easy  to  understand,  of  leading  to  re- 
latively simple  analysis  and  computation,  and  of  giving  a  scattering  cross  section 
which  agrees  with  that  given  by  the  spheroidal  function  solution  at  intermediate 
frequencies. 

It  is  desirable  to  have  a  similarly  simple  approach  at  low  frequencies  since 
all  the  methods  described  require  complicated  mathematics  and  considerable  computa- 
tion.    A  classical  nethod,  due  to  Rayleigh,  based  on  a  static  approxiriation  requires 
only  simple  analysis,  but,  unfortunately,  has  only  a  very  restricted  range  of  appli- 
cability.    The  work  that  follows  is  an  attempt  to  simplify  as  many  aspects  as  possi- 
ble in  the  low-frequency  problem.     A  comparatively  simple  fiiysical  picture  is  employ- 
edj  the  ensuing  analysis  is  complicated  but  enables  one  to  reach  a  stage  where  the 
computation  is  simple.     The  matherre  tical  technique  has  been  developed  only  for  the 
circular  disc  but  it  is  hoped  that  it  can  be   generalized  to  more  general  obstacles. 

The  physical  picture  employed  is  a  generalization  of  Rayleigh's  static  picture. 
In  describing  it,  it  is  convenient  to  use  the  terminology  of  sound  waves.     At  low 
frequencies  the  field  is  essentially  static  in  any  domain  whose  linear  dimensions 
are  much  smaller  than  a  wavelength  and  there  is  very  little  radiation  from  the  disc. 
If,   therefore,  we  solve  the  problem  when  there  is  no  radiation  we   should  obtain  a 
very  good  approximation  to  the  pressure   distribution  on  the  disc.     We  can  then  con- 

struct   a  radiation  field  which  has  the   same  pressure   distribution; this  should  be 

Clemmow  discusses  electromagnetic  phenomena,   but  the  fundamental  principles  are 
applicable  to  scalar  problems. 
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an  excellent  low-frequency  approximation  to  the  field.     Clear iy  we  could  devise 
an  iteration  scheme  based  on  these  ideas  which  would  improvs  the  accuracy  and  ex- 
tend the  frequency  range  of  applicability.     The  difficulties  involved  in  the   dis- 
cussion of  the  convergence  of  an  iteration  scheme  are  well  known  and  it  is  there- 
fore desirable  that  the  mathematical  formulation  shall  be   such  that  these  difficul- 
ties either  do  not  occur  or  can  be  overcome. 

It  must  be  realized  that  the  problem  in  which  there  is  no  radiation  is  not  the 
static  one;  it  is  the  problem  in  which  the  waves  at  infinity  are  standing  waves. 
The  approximation  of  standing  wa-ves  is  much  less  drastic  than  the  static  approxima- 
tion because  the  field  still  satisfies  the  harmonic  wave  equation  although  It  is 
non-radiating. 

Accordingly,  in  the  mathematical  formulation,  we  split  the  scattered  field 
into  a  non- radiating  field  which  becomes  the  static  field  in  the  limit  of  zero 
frequency,  and  the  remainder  of  the  field.     This  remainder  is  then  added  to  the 
incident  field  and  the  svim  is  treated  as  the  incident  field  in  the  non-radiating 
problem.    This  problem  is  then  solved  exactly  by  means  of  integral  equations  and 
leads  to  a  Fredholm  integral  equation  governing  the  original  problem.     In  this  way 
the  question  of  the  convergence  of  the  earlier  iteration  scheme  is  avoided.     The 
iteration  is  obviously  equivalent  to  solving  the  Fredhola  equation  by  a  Neumann 
expansion,  so  that,  if  we  wish  to  use  it,  we  can  quote  well-known  theorems  to  segr 
when  it  will  converge. 

The  mathematical  difficulty  mentioned  earlier  arises  in  the  exact  solution  of 
the  non-radiating  problem  and  the  method  can  be  generalized  to  other  shapes  only  if 
this  solution  can  be  obtained  in  general.     It  should  be  remarked  that  the  method  can 
be  used  for  the  diffraction  of  electromagnetic  waves  by  a  circular  disc;  this  will 
be  discussed  in  a  future  paper. 

The  Fredholm  integral  equation  has  many  advantages.     The  power  series  solution 
can  be  obtained  from  it  very  rapidly.     It  is  easy  to  solve  approximately;  one  ob- 
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tains  a  pressure  distribution  vhich  can  be  computed  easily  from  well  tabulated 
functions.  Further,  it  is  possible  to  estimate  the  error  made  in  solving  approxi- 
mately a  Fredholm  integral  equation,  so  that  one  can  state,  for  example,  that  our 
value  for  the  scattering  coefficient  is  correct  to  so  many  percent.  Also  it  is 
possible  to  deduce  an  upper  bound  for  the  scattering  coefficient.  A  Fredholm 
integral  equation  has  been  derived  by  Magnus"-  -I,  but  it  is  much  more  complicated 
than  the  one  wa  obtain. 

The  following  anailysis  is  restricted  to  the  case  when  the  pressure  distribution 
on  the  disc  depends  only  on  the  radius.  The  resulting  integral  equation  is  ^'ormu- 
lated  in  Section  1.  In  Section  2  the  angular  integration  is  performed,  and  in 
Section  3  the  Fredholm  integral  is  derived.  The  problem  of  the  incident  plane 
wave  is  discussed  in  Section  U  and  an  approximate  solution,  agreeing  well  with 
the  spheroidal  fiinction  solution,  is  obtained  in  Section  5.  Section  6  is  devoted 
to  the  problem  when  the  incident  wave  is  supplied  by  a  point  source. 

In  Section  7  w©  give  a  general  theory,  applicable  to  any  obstacle,  which  shows 
how  an  upper  bound  for  the  scattering  coefficient  may  be  obtained.  Section  8  con- 
tains an  estimate  of  the  error  introduced  by  our  method  of  approximation. 

The  derivation  of  a  number  of  results  has  been  carried  out  in  appendices. 
References  to  equations  in  appendices  are  written  for  example,  as  (Bl),  meaning 
equation  (1)  of  Appendix  B» 

1,   The  integral  equation 

We  consider  the  problem  of  the  diffraction  of  a  harmonic  sound  wave  by  a  rigid 
circular  disc  of  radius  a.  The  boundary  condition  is  that  the  normal  derivative  of 
the  field  must  vanish  on  the  disc.  The  case  when  the  field  vanishes  on  the  disc 
could  also  be  treated  by  the  method  we  shall  give  but  will  not  be  considered  here. 
Results  for  the  circular  aperture  in  a  plane  screen  nay  be  deduced  from  those  for 
the  circular  disc  by  Babinet's  principle. 
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Let  p,  0  and  z  be  cylindrical  polar  coordinates^  and  normalize  all  lengths  by 
division  by     a     so  that  the  disc  occupies  0<p<l,  z-0.     Let  the  incident  field 
be  independent  of  (2  and  be  represented  by  u  (p,z)  at  (p,5l,z),  the  time  dependence 
e  ^    being  understood.     The  total  field  u(r),  at  the  point  whose  position  vector 
with  respect  to  the  origin  is  r,   is  then  independent  of  (jS,     Such  a  field  occurs 
for  a  normally  incident  plane  wave  and  when  the  incident  wave  is  due  to  a  source 
on  p  ■  0»     The  problem  of  the  unsymmetrical  incident  field  will  be  considered  in 
a  separate  work. 

The  total  field  at  the  point  r.   can  be  represented  by 

u(r^)     .    ujp^,z^)  .  ^    I  g(p)   j"^T(r,r^)]        dS 

-ia|r  -r|  / 
where  S  is  the  area  of  the  unit  circle,  p  is  a  point  of  S,  T  *  e        /|r--r|  , 

and  a  is  the  (dimensionless)  product  of  the  wave  number  and  a.  To  ensure  a  unique 

solution  (Jones'-  -' )  we  have  to  impose  the  edge  condition  that  g  »  0  on  the  boundaiy 

of  S. 

The  boundary  condition  that  has  to  be   satisfied  on  S  is  9u/3z,   -  0.     To  avoid 

non-integrable  singularities  in  the  integral  we  modify  the  expression  for  3u/9z, , 

We  have 


9u 

az 


'-  ■  sir  W^)  ♦  sir  U'a)  rsfc^'^'^-i']  „ "' 

11  /g  '-     1  -Jz-O 

•  st^^o'"!"!'  *  Er|  ^(£'  fe  %^  *  "j  ''^'^i>  "= 

2         2  2 

since  3Y/az  ■  -aY/3z-   and  (V     +  a  )  T  »  0^ where  V     is  the  three-dimensional  Laplacian. 

Since  g  ■  0  on  the  boundary  of  S  the  divergence  theorem  gives 
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~  ■  ^%^^1'^1^  *  117  I  V  g(p)T(p,rjL)  -  grad  g  .  grad  T(p,iVj^)LdS 


32 


where  the  operators  grad  and  div  are  independent  of  z,  the  suffix  1  indicating 
operations  with  respect  to  p  »  There  is  now  no  difficulty  in  taking  the  limit 
as  z^  ->  0  and  we  obtain  the  following  equation  to  determine  g 


(1) 


■  al"  %^Pl»°^  "  B  '^^''l  j  '^^e.'Pi^  grad  g  dS  +  J^  j  g(p)T(p,p^)dS, 


where  p,  is  a  point  of  S. 

Since  g  is  independent  of  {f,  only  T  is  involved  in  the  integration  :7ith 
respect  to  (jS,     We  proceed  to  the  consideration  of  this  in  the  next  section. 

2»       The  integration  with  respect  to  gf 

In  the  first  place  we  shall  consider  only  the  real  part  of  f.     Then,  it  is 
necessary  to  evaluate  integrals  of  the  type 


/^"  cos  r  \jl  +  b^-2b  coe(0-9(  )[  ■'■^^ 

I  -     cos(gr-?)    \ ^ d9(. 


jl  +  b^-2b  cos((2-5(j^;  I  •'■/^ 


where  p,  y  ^^^  b   (<  l)   are   arbitrary  real  constants,   since 


|p  -  P^l  -    ^  P^  *  Pi  -  ^PPi  cosCgr-gr^)  ^  ^^^ 


Now 
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U  co8((7  -p)       cos  2gf     Li 2Li_-_  d0 

•^        2  (i+b  -2b  co52(2)-^^ 


I 

■  X  ■ 

'o 

2  cos(gf^-p)    f    1 
o 


on  putting  tan  9  =  cos  Sg-b   *     f^®"°® 


cose(l-b^sin^©)-'-/^+bsin^©|  cosY|(i-b^8in^9)-'/^-bcos©| 


(1-h^  sin^©)^/^ 


b/  "'"  ■*■  1  1  ' 

r5^     T  -  **  ,.r.^(ri    >^^    f      x^cos  yd-x^)^  cos  yCb^-x^)^        •       /n     2,Z,        .^2     2,Z 
(2)    I  -  ^  cos(gr^-b)  J  j         (i.x^)V^  (b^.^^j^/^ —      *  '^"  ""  "" 

Similarly 

z^"  cos  Yfj-+b^-2b  cos^^-grjU/^  /^  ,,      2.1/2  ,^2     2,1/2 

-To^      I  L ^1 J  j/K  _  1    I     cos  y(1-x  )   '    cos  Y(b  -X  )   '        , 

Similar  resvilts  may  be  obtained  for  the  integration  of  the  imaginary  part  of  Y 
with  respect  to  p  but,  for  the  form  in  which  we  want  them,  it  is  quicker  to     proceed 
in  a  different  fashion.     We  have 


dx. 


where  J    is  the  Bessel  function  of  the  first  kind  and  v-th  order.     Cn  using  Sonine's 
first  finite  integral  (Watson ^    J)  the  integral  on  the  right  may  be  written 


2n       n/2 
a 

o        o 


I         (sin  ©  J^(a|p-p^|sin  ©)  d©  d(jf 


n/2 

1" 


•  2na   I     sin  ©  J  (opsin©)  J  (ap,sin©)d© 


after  a  use  of  the  addition  theorem  for  Bessel  functions.     Hence 
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sinalp-p^l  2    /     ^  sina|p-Pj^| 

grae   -   -" 

it/2 


S  S 


j  pg'(p)  2na     I       sin  6  J-j^(apsin9)   J^(ap-j^sine)d&dp 


1  n/2 


+  2na     j   Pg(p)    I   sin©  J^(apsine)  J^(ap, sine)d9dp 


'o  'o 

1  n/2 


■  2iia         Pg(p)         sine  cos  6  J  (opsin©) J  (cp, sin©) d©dp 


(U)  2Tia^       d 


o  o 

1  n/2 


~-    ^       P-j^Pg(p)         cos  ©  jQ(apsin©).Jj^(sp-j^sin&)d©dp   . 

0  0 

Now,  on  putting  \  ■  1,  ^  -  j-  and  v  ■  «■  in  the  formula  (C.l)  given  in  Appendix 
C>  we  have 

"/2                         Z    ^    f/„2  2x1/2]   -    f,  2  2sl/2) 
(5)     J  (zsin©)JT(Zsin©)coE  ©d©  -  S-   x    ^^  L  -   2_   ^^^  S  i  ;i ^  dx. 


;ince 


) 


the  use  of  (2),   (3),   (U),   and  (5)   in  (1)  leads,  after  some    straightforward  mani- 
pulation,  to  the  equation  ,                               ^ 

1  inin(p,p, )   ^        /       p     ?  Ti           {  ,   2     2^) 

,        H      /     .  {               ^  t^cos  a(p^^-t^)'^Uos  aCp'^-t^)^ 

3z^     0     1             np^  dp^  j  j                             (p2.t2)l/^  (p^-t^)^/^ 


dtdp 


^1^)°'^*'']    lOT7r-3r  77:^5^175— 

0  0  1 


1 

f 

dt 
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or 


dz.      o     1 


"Pj     dpj 


(6) 


(    t^g'(p)cos'fa(p^-t')^4    .„ 


d       »lnW:iVf.l; 

3r    — -S — ^  T  ,6 'op 


4    f  ^4.    c-^    cl       sin{a(p 


The  solution  of  (6)  is  considered  in  the  next  section. 


dt  . 


3. 


The  solution  of  the  transformed  equation 
The  equation  (6)  is  of  the  form 

P 


fCp^) 


dt 


The  solution  of  an  equation  of  this  type  is  given  in  Appendix  A  (see  (A.U)). 
Identifying  h(t)  with  its  appropriate  counterpart  in  (6)  wb  have 


or 


(7) 


(a(p2.x2)V2} 


/     ,         cos|a(p  -X  }  '  j  /  sinJa(p  -x  )   '   ^ 


in(a(p2.x2)V2J 


-  -2  f  ^  u 


(t) 


t  cosh  I  aCx^-t^)-"-/^ 


TT^t"^?^ 


dt 


«  /   a        *  cosh  ^  a(x  -t  )  M 


There  ie  a  number  of  vaye   of  rewriting  (7).  We  will  consider  the  method  which 
is  applicable  to  low  frequencies  or  small  values  of  a.  The  first  tern  of  the  left- 


-  10  - 


hand  side  was  derived  from  the  part  of  the  original  integral  equation  which 
corresponds  to  standing  waves  at  infinity  and  which  goes  over  to  the   static 
limit  as  a  ->  0.     Following  the  principle  described  in  the  introduction  we 
regard  (7)  as  an  equation  in  which  the  first  term  on  the  left-hand  side  is 
equal  to  a  known  function.     The  solution  of  such  an  equation  with  the  condi- 
tion g(l)  -  0  is  given  in  Appendix  A  (A, 5)   and  applying  this  solution  to  (7) 
we  obtain 


g(x) 


2    /    t  cosh|a(t  -X   )   '  I  2    / 


p  c0Eh{a(t^-p^)^/^ 


°^'^  (t^.p2)V2 


dp 


(8) 


-  i 


.1           ,       sin  L(p2-t^)^4 
^'P'  dp     ,\  J,l/l ^ 


(p5:??75 — '' 


dt 


Thifl  equation  holds  for  0  <  x  <  1,     Now 


J. 

1 


[ait'.^)'^^] 


t  cosh  <  a(t  -X  ) 


(t5-.^)V2 


g(p) 


dp  dt 


.2il/2 


on  inverting  the  order  of  integration  and  making  obvious  changes  in  the  integrating 
variables.     Also 

,y  t  slnL(tV)'/'}  cosh  Uj^-t^)'-"'] 

7     sinh{a(p^-t^)l/^  .  a(ir^-t^)^''^}  *     slnh  {a(p^-t^)^^^  a(y'-t^)^^^} 


I 


1 
5" 


tdt 


-  il  - 


Ap  -y  ) 


(9) 


2  < 


sinh  at  ,. 
7 at 


p+y 


sinh  at 


dt 


P-y 


Hence  (8)  becomes 


(10) 
where 

(11) 

and 
(12) 


G(y)  -  F(y)  ♦ 


i    f    q(3j)  J    sinh  a(x-y)       sinh  a(x4y)    I 
«  J  I  x-y  xTjF  r 


dx  (0  <  y  <  1), 


G(y) 


(i.y^)^/^ 


.  ?-  coshJa(t^-l+y^)^/2l 

{l.y2)l/2    (t2.1+y2)^/2 


14" 


cosh/a(t^-p^)''-^l 


°''"    (t^-p2)Vi 


dtdp  • 


It  should  be  remarked  that  the  method  of  deriving  (10)  shows  that  the  unique 
solution  of  (1)  is  a  solution  of  (ID).  Since,  however,  the  eigenvalues  of  (10) 
are  real, the  theory  of  the  Fredholm  integral  equation  implies  that  (10)  possesses 
a  unique  continuous  solution.  The  solution  of  (10)  can,  therefore,  be  identified 
with  that  of  (l). 

We  shall  nov7  consider  equation  (iO)  for  two  particular  incident  fields, 

U.   The  incident  plane  wave 

Suppose  now  that  the  incident  field  is  a  plane  wave  so  that 


u^(p»z)  -  e 


-laz 


Then 


g|-  u^(p,0)  -  -ia 


and  we  obtain 
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.2_-,   2,1/2: 


,                      cosh  4  a(t  -1+y  )    , 
^^^^"T         2  1/2   b,  ^1/2  ^  ^^  ^^  ^^ 


-^1 


l+y 

sinh  ax  , 
—IT-  ^ 

1-y 


in  the  same  way  as  was  used  in  deriving  (9).  Hence  (10)  becomes,  for  an  incident 
plane  wave. 


1+y  1  \ 

(13)       G(y)  -  -  ^  f      ^^^^^-^  dx^ll  Gi^)   hlnh  ajx.y)   ,  sinh  a(x+y)K^  (0  <  y  <  i) 
n    ;             X                   n   I             1        x-y  x+y  f  -  •'   - 

l-V  ft  I  1 


Observe  that  the  right-hand  side  of  (±3)  vanishes  when  y  •  0,  as  it  should. 
In  fact  the  right-hand  side  is  an  odd  function  of  y.     Therefore  we  can  define 
0(-y)  when  y  is  positive  to  be  -G(y)j  equation   (13)  then  becomes 


(lU)  G(y)  -  -  ^  £inh_ax  ^  ^  i 

n    J  X  n 

1-y  -1 


1+y  1 

(    G(x)£2£\4i^dx     (-l<y<l). 


This  Fredholm  integral  equation  has  to  be  solved  to  complete  the  solution  of 
our  problem.  There  are  several  methods  available  and  the  application  of  some  of 
them  will  now  be  discussed. 

Firstly,  an  exact  solution  of  (lU)  may  be  written  down  by  means  of  the  Fredholjn 
theory,  but  the  form  of  the  kernel  makes  the  calculation  of  the  resolvent  impractical. 

Secondly,  we  could  find  a  solution  by  expanding  in  terms  of  the  eigenf unctions 
of  the  kernel.  The  eigenfunctions  are,  however,  spheroidal  functions  (Meixner  and 
Schaf ke "-  -'),  and  so  this  method  appears  to  have  no  advantage  over  separating  the 
harmonic  wave  equation  in  spheroidal  coordinates. 

Thirdly,  we  could  attempt  to  form  a  solution  by  means  of  iteration.  The 

1  ATT 

simplest  way  of  doing  this  is  to  take  -  —  I —  dx  as  the  first  approximation 


X 

1-y 
and  substitute  successively  in  the  integral.  This  is  an  excellent  method  when  a  is 
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small  because  the  kernel  starts  at  the  term  involving  a     since  G(x)  is  odd.     Hence, 
since  the  first  approximation  starts  at  a,  the  first  approxination  to  G(y)  is  correct 
to  0(a),  the   second  to  0(a  )  and  so  on.     It  is  easily  verified  that  the  expansion 
in  powers  of  a  is  the  same  as  that  obtained  by  Bouwkamp  by  substituting  a  power 
series  expansion  for  g  in  (1),     However,   the  iterated  solution  of  (lU)  provides 
the  coefficients  in  the  expansion  much  more  rapidly  than  his  method.     Two  iterations 
can  be  carried  out  very  quickly.     Since  the  resulting  pressure  distribution  contains 
more  terms  than  given  by  Bouwkamp  it  is  worthwhile  to  reproduce  it.     It  is 

r(^r)     -     -  !ii22.   Ji  ^  <^^       2ia^  +   °^        13ia^  ,    /I  U    \  6       323iJ 


Sin 


/I  36     V  8 


2ia3-3 


£-Z- Ji  +  3a     .  2ia^  ^  a'*  .  19ia^   .    /  1  h    \6 


2ia^y^  J,   ^  5a2       iOia^  ^  5a'*  ^ 


£ia!z!     i,7a\  2iaV 

I76I4DH  r  *  iF"  *    •••    (      1632960n 


For  comparison  with  Bouwkamp 's  formula  it  must  be  remembered  that  his  distribution 
is  for  the  complementary  aperture  problem  and  is  therefore  -  «■  times  ours,  also  he 
uses  the  time  dependence  e         , 

It  is  known  that  the  iterated  solution  described  fails  to  converge  when  1/n  ex- 
ceeds the  first  eigenvalue,  i.e.,  it  fails  to  converge  when  a  exceeds  some  constant. 
This  probably  means  that  the  power  series  expansion  also  diverges  when  a  exceeds 
this  value  but  it  may  mean  only  that  the  expansion  has  been  rearranged  in  a  diver- 
gent series  of  groups  of  terms.     The  iterated  solution  converges  if 


-  lU  - 

(X5)  \    [     f     £i2*L'_jiii:lI   dx  dy    <    1, 

which  is  certainly  true  if  sinh  2a  <  n.     It  does  not,  of  course,  follow  that  the 
power  series  expansion  converges  for  this  range  of  a. 

More  complicated  methods  of  iteration  are  available  which  converge  everywhere 
(Bruckner L    -I)  but  no  attempt  has  been  made  to  use  then  from  equation   (lU).     Instead 
we  have   solved   (lU)   approximately  in  a    straightforward  manner  which  leads  to 
tabulated  functions,  enables  G  tc   be  computed  for  moderate  values  of  a  fairly 
rapidly^  and  provides  an  estimate  of  the  error. 

5.       Solution  by  approximation  to  the  kernel 

Since    (8inhaz)/z  is  an  integral  function  it  can  be  expanded  in  a  Taylor  series 
about  any  point.     We  deduce  at  once  that 

1+y  ^     „     2n  2n+l 

oj    /  •    u  00     a  a    y 

„/    \  2i    I         sinh  ax    .         r—      n       •' 

G^y)     -     -T-j         -ir-^*r  (2n.l)l     • 

i-y  "-^ 


The  substitution  of  this  representation  in  (iU)  leads  to  the  equations 

I       2m  _ 

nia  a  oo 

(^6)  *     r^iTT    "  -  ^iV     ^     i:    ^mn  ^n  ^^  "  °'  ^>    '"^ 


nK) 


to  determine  the  constants  a  .  Here 

m 


2n  /       o     1      j2m+l       .    , 

,  a  I       2n-t-l     d  sinh  ax    . 

^mn  ■  (2m+l)l(2n+l)i     I     ^  ^TliJTT    — 3E ^ 


and 


\     '     (2m+i)l 


1     l>x  2jn+i 

''     f        sinh  at     d  sinh  ox      ..    ^ 

t  ,^2m+l  x 
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It  is  shown  in  Appendix  B  that 


*inn"  ^m"   (2m+l)  U2n+i;i(n-m; 


n      2m+2E 
a 


m        2n+2s 


''"*'"  St^  ■  *'"*'"  &  ^^^  r"" 


n-l     2in+2s+l 
a 


m-1     2n+2E+l 


E=0 


s=0 


(m  /n) 


and 


2in 


2     1 

mm     ( 2ra+i ) I 


^i"^  °  r     (f^^   Y(2iB+2)   -  T(2s+1)  j,  -  ^  I   Fi(a)  -  Ei(-a) 


m-l       2s+l 
a 


-  cosh  a  J2    ^JiTXTTS  ^(2in+2)  -  Y(2s+2) 


where  T(bj)  =  1  +  ^  +  ...  ♦~j^-  Y»  Y(l)  -  -  y»  Y  "  0.5772   ...  and 

Ei(x)     ' 


^       e* 

l-dt. 


^-00 

The  scattering  cross  section  of  an  obstacle  with  an  incident  plane  wave  is  -(Uir/kJv/A 
where  A  is  the  amplitude  of  the  scattered  wave  in  the  direction  of  propagation  of  the 
incident  wave  and  k  is  the  wave  number.     Thus  the  scattering  cross  section  of  the  disc 

^=  1 

-2Ba    ^     I     gpdp  . 

Hence  the  scattering  coefficient  (ths  cross  section  divided  by  the  area  of  the  disc) 
is 


a  a 

n 


2n 


00 

-2^  ^  (2n*i)l(i!n+3) 


The  complete  solution  of  (16)  is  scarcely  practicable  but  an  approximate  solu- 
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tion,  in  >diich  the  infinite  series  in  G  is  truncated  to  form  a  finite  series,  can 
be  provided  with  reasonable  computation  when  the  tnincated  series  does  not  involve 
too  many  terms.  A  method  for  calculating  the  error  produced  hj   this  truncation  is 
given  in  Section  8« 

There  is  no  possibility  that  the  determinant  of  the  coefficients  vanishes  either 
for  the  complete  or  truncated  series.  If  the  deterndnant  did  vanish^ i/n  would  be  an 
eigenvalue  of  the  kernel.  This  is  impossible  since  the  eigenvalues  of  both  the  com- 
plete and  approximate  kernels  are  real. 

In  particular,  on  taking  only  the  first  term  of  the  series  in  (16),  ve  obtain 

2iK 


00 


^^^  \o   ■  ^  ^^^^  °-  '   ^^^°^  *   Ei(-a),  and  K^  -  j-Je'^EiCUa)  +  e°Ei(-Ua)-2cosha(r+lnaa) 


The  scattering  coefficient  is  then 
UK 


0 

3(A^     +  /) 

^00 

We   show     in  Table  1     some  results  obtained  by  using  the  approximation  in  which 
it  is  assumed  that  only  a     is  non-zero.     The  approximate   scattering  coefficient  is 
correct  to  within  about.  10 /j  for  o.  <  2  but  thereafter  diverges  rapidly. 

Table  1 


a 

coefficient 

0.5 

0.00200 

1.0 

0.0387 

1.5 

0.253 

2.0 

O.BBli 

2.5 

1.617 

3.0 

2.039 

One-half  the  scattering  coefficient  of  a  circular  disc  when  a     -  0  for  n  >  1. 

n  — 
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The  scattering  coefficient  obtained  by  assuming  that  only  a     and  a,    are  non- 
zero is  given  in  Table  2,     The  values  obtained  from  the  approximation  were   calculated 
correct  to  U  significant  figures.     Also  shown  are   the  correct  values  of  Bouwkamp, 
It  can  be  seen  that  the   approximate  values  agree  with  Bouwkamp's  to  within  U/o  for 
a  <  2,5  (within  ^   /o  for  a  <  2),     (For  an  estimate  of  the  error,   independent  of 
Bouwkamp's  results,   see  Section  8  and  Table   6.)     The  value  at  a  «  3  is  13/5  too 
high}  this  can  be  improved  by  taking  more  terms  in  the   series.     Thus  if  we  take 
a_  non-zero  we  find  that  the  scattering  coefficient  is  2  x  1.158,  which  is  correct 
to  within  1  »•  /o» 

Table  Z 


a 

Correct 

Approximation 

a 

Correct 

Approximation 

0.1 

0.00000300 

0.00000301 

1.6 

0.38U 

0.3819 

0.2 

0.000CU86 

O.OOOOU865 

1.7 

0.505 

0.5019 

0.3 

0.000250 

0.0002503 

1.8 

0.6U2 

0.6378 

0.U 

0.000809 

0.0008089 

1.9 

0.783 

0.7812 

0.5 

0.00203 

0.002032 

2.0 

0.922 

0.9210 

0.6 

0.00U36 

O.OOU363 

2.1 

1.0U2 

1.0U6 

0.7 

0.008U3 

0.008U22 

2.2 

1.135 

1.1U8 

0.8 

0.0151 

0.01506 

2.3 

1.200 

1.22li 

0.9 

0.0255 

0 .025U1 

2.U 

i.235 

1.27U 

1.0 

O.OiOl 

o.ouioo 

2.5 

1.2U8 

1.3016 

i.l 

0.06U 

0.06380 

2.6 

1.2/ili 

1.313 

1.2 

0.097 

0.09626 

2.7 

1.227 

1.312 

1.3 

0.1ii2 

0.1U13 

2.8 

1.202 

1.303 

1.U 

0.203 

0.2022 

2.9 

1.173 

1.290 

1.5 

0.283 

0.28l5 

3.0 

l.lli2 

1.272 
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The  column  headed   'Correct'  is  reproduced  from  Bouwkamp's  review  and  gives 
the  transmission  coefficient  of  an  aperture.     The  column  headed    'Approximation' 
is  one^half  of  the  scattering  coefficient  of  a  disc  determined  by  assuming  that 
only  a     and  a^  are  non-zero.     The  value  corresponding  to  a  ■  3  when  a  ,  a.,   and 
a^  are  non-zero  is  1,158. 

One  point  that  should  be  observed  is  that  the  approximation  reproduces  the 
maximum  which  occurs  in  the  correct  values  although  the  magnitude  is  some  5/o  in 
error. 

For  small  values  of  a  the  approximations  (B13)   -  (B17)  inclusive  can  be 

used.     The  corresponding  approximations  to  a     and  a,   are 

(18)  /^(a^)   ^    .     -^    (^  .  ^),    ^(a^)  ^    --^  .  —^  , 


(19) 

/^(V 

a    - 

oa 

i 

The  1 

scattering  coefficient  is  i 

twice 

8.'' 

(1* 

8      2 
25"° 

approximately. 

6.   The  symmetrical  point  source 

We  now  consider  what  happens  when  the  incident  field  is  due  to  a  point  source 
on  the  axis  of  the  disc.  Suppose  that  the  point  source  is  at  (0,0,-8  )  where  z  >  Oj 


then 


^^^*^1^  '  (2   .„     _  ,2)  1/2 
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Since 


\  3z,  /  \dz 


,  we  have 


Z--0  2,  "0 


1  11 

gi^     ^  cosh  [a(t^-l^y^)^j  .*     cosh  [a(t%^)^J    exp  i-iaCp^+z^jH 

anr  J        ,.    ,  ,5  ..2  ,      2.1/2  P       ,^2     2,1/2 ,   2     2,1/2 ^ 

(1-y  )^    ^  ^  ^^  ^  o    (t  -p  )  '      (p  +Zq)  ' 


or 


1  (2    2  2') 

,   -iaz   /        cosh^a(t  -1+y  )  \   tcosh  at  +  iz  sinh  at 

By  the  same  method  that  (9)  was  derived  this  may  be  converted  into 

1+y     2    2 

2  ~la2^  /  '   2(t  +l-y  )ccsh  at  ♦  Uitz  sinh  at 


(21)  F(y)  -  -  1  e         " ^ dt. 

l.y  (tSl.y^)^TV^ 

The  integral  equation  for  G  may  now  be  written  down.     The  solution  of  this 
equation  can  be  used  to  obtain  the  scattering  coefficient,    (which  is  defined  as 
the  ratio  between  the  average  rate  at  which  energy  is  scattered  and  the  average 
rate  at  which  energy  is  incident  on  the  disc);  the  result  is 
(l,,2)l/2  1  ia(p^z?)^/2 


(22) 

a) (1*2!) 


(,-,     2,1/2 ]    ^  \     ^  di ,  .2     2,1/2   ■  P  **P' 

^(l*z^)       -z^j       ;  0       (p  +z^) 


(For  a  general  proof  see  Jones'-    J.) 

Vfe  shall  n 
(20)   shows  that 


Vfe  shall  now  restrict  our  attention  to  the  case  when  z_  is  large.     When  z     »  1, 

0  o  ' 
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-ia» 


F(y) 


[ 


where 


F^Cy) 


l*y 

2i      f  sinh  at    ,^ 

T    J  — t — ^** 

1-y 


^2^y^  "  "  i?^   ^^"^  a(l+y)-sinh  a(l-y)*a(l+y)co8h  a(l-y)-a(l-y)coeh  a(l>y) 


1+y 

1-y 


*  a^d-y^)    I         sinh  at  ^^ 


and 


F3(y) 


k[ 


■L+y)  ^  J  -  ?  a(l-y)^J'  cosh  a(l+y). 


\*  ^l-y)^leinh  a(l*y) 


1      1 


-^^■y^     ia'2  «(l*y)7  c°^h  a(l-y)   *  ^   ^  *  ^l^y)^).  sinh  a(l-y) 

i+y 


+2(l-y^)^l  +  i  a^(i-y^) 


sinh  at 


dt 


L-y 


The  reason  for  including  the  term  involving  F^  will  be  given  below. 


Let 


-iaz 


G(y) 


0_^(y)^^G2(y).i^G3(y).    0(^^)1 


Then,  since  F(y)  is  odd,  the  equation  (10)  may  be  transformed  into  the  form  (lU)j 

hence,  taking  corresponding  powers  of  —  together,  we  obtain 

^o 


(23) 


X 

G^(y)  .  F^(y)*^  I  G^(x)  ^inh_a(x-y)  ^^     (.1  <  y  <  1) 

-1 


with  n  -  1,  2  or  3. 
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In  view  of  the  expression  for  F,,  G,   satisfies  (lU),  i.e.,  it  is  the  solution 
for  the  incident  plane  wave.    With  regard  to  Go  and  G^,  we  put 


and 


,      2n  2n+l 
00      bay 

GgCy)    -    FgCy)  *    E       (2d.1)  t 

n-o       ^  ' 


^  „     2n  2n+l 

CO  cay 

G3(y)     -     F  (y)  .    Z        "(W)i         • 

n-o 


The  equations  for  b     and  c  ,  analogous  to  ( 16),  are 

,,      2m 
tiib  a  00 


^^^^  -ronryr  "  \  *  ^  V\  (m-o,  i,  ...) 

and 


n»o 


(25) 


where 


2in 

nic  a  00 

m 


T^TTTT    -    \     *    C      V^n  (m-0,  1,   ...) 

^  n»o 


J  1 (      ^  f.  d^^      sinh  ax     . 

•  ^-1  °^ 


and 


jj.  ox 

The  expression  for  the  scattering  coefficient  may  also  be  approximated  when 
z    »  1.     In  order  to  remove  an  exponential  from  the  integration  we  shall  make  the 
additional  assumption  that  a/2z    «  1.     Then  the  scattering  coefficient  is 
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-2/e 


O^y  dy 


--^i 


az 


ia  G.  -  <Ii  +  ^a^(l-y2)y  G^ 


o      — 


y  dy 


(26) 


2    J 
az 


o     o 


iaG,  +  p  ioG,    - 
3       U         1 


|1  +  ^  a2(l-y2)l  G2.ia  ]^  (l-y^)  ♦  ^  (1-y^)^  G^ 


y  dy 


+  O(^)  . 

z 

o 


The  first  term  is  the   scattering  coefficient  for  a  normally  incident  plane  wave. 
We   shall  now  prove  that  the   second  term  is  identically  zero,   so  that  there  is  no 
correction  to  the  plane  wave   scattering  coefficisnt  if  the  term  F-  in  F  is  neg*- 
lected.     In  fact,   it  will  be   shown  that  the   scattering  coefficient  contains  only 
even  powers  of  l/z   . 


The  integral  equation  (1)   can  be  written  as 


(27) 


r3u    -, 


L 


TG 


where  T  is  a  linear  operator.  In  (27)  the  brackets  indicate  the  value  at  z,  »  0 

2 

and  ve   imagine  that  9u  /dz  is  expressed  in  terms  of  1-p  .  We  msy  then  write  the 

integral  in  the  scattering  coefficient  (22)  as 


r3u  -t 


y  dy 


(say). 


where  the   star  indicates  a  complex  conjugate.     Clearly     [  G,       -rr-       )"  I      tt"   L 


Suppose  now  that,  when  z     »  1  and  a/2z     «  1, 
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[si^]-    E.'X  *0U/.«) 


-iaz. 


and 


N-1 


£i  v^o  ^  0(1/ 


where  v  and  G  are  independent  of  z,.  Then,  from  (27), 


n      n 


and  the  integral  I.  in  the  scattering  coefficient  is  given  by 


Now  Vp  is  purely  real  and  v^  ,  is  purely  imaginary.  Hence 


n-1     ^  n-1 


7~  (G  ,v   )  -  r  6   (G  ,v  „) 
■^  ^  p*  n-p     ^  n-p  p'  n-p 


p=l 


p=l 


P  P  n-p 


■ 

g  Vp(°p.^«n-p> 

m 

n-1 

y~    5       (TG  ,G       ) 
L^     n-p^     p»  n-p^ 

m 

r  Vp('p'Vp^ 

n-1 


)~  6  6   (G  ,v   ) 
<--,  p  n-p^  p*  n-p' 


p-1 


where  6  ■  ±  1  according  as  p  is  even  or  odd.  Since  5_  6 

p  P  n— T 


P  n-p 


±   1  according  as 


n  is  even  or  odd  we  see  that 


1  2n-l      . 

ii  -  E    4  E  %-  '*„-p>  • 

n»l  z    P"0 
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Further,  the  factor  of  I_  in  (22)  involves  only  even  powers  of  z  and  hence  the 
scattering  coefficient  contains  only  even  powers  of  1/z  .  This  proves  our  earlier 
assertion* 

The  third  term  in  (26),  therefore,  is  the  first  which  provides  a  correction 
to  the  plane  wave  scattering  coefficient.  This  is  the  reason  why  it  is  necessary 
to  include  the  term  F.  in  the  expansion  of  F  in  inverse  powers  of  i   .  The  scatter- 
ing coefficient  for  the  point  source  when  z     »  I   and  a/2z  «  1  is  the  plane  wave 
scattering  coefficient  together  with 


t        n-o 
o 


2n 
g 

(2n+l)l(2n+3} 


2   b      ,, 

1  c  -  (1+  ^)   ~  ♦  i  Vf 
n     2n+5  a     \h 


(2n+5)(2n-.7) 


n 


*  0(i/i!i). 


There  is,  of  course,  no  contribution  from  the  terms  involving  Fp  and  F^  in  the  ex- 
pressions for  Gp  and  G-  since  these  are  derived  from  fields  which  do  not  radiate. 
Instead  of  solving  (2U)  and  (25)  we  can  find  G^  by  means  of  iteration  when 


a  satisfies   (15)    •     The  power  series  expansion  for  Gp  after  two  iterations  i« 

8ln 


r     .       2y  Jp  ^     2      h±  J  ^  cT      U6ia^  .    /  1  8    \     ^ 


13981    7  ^ 


[m 


1H9U0 


96     \     8  ^ 


2025n 


^  2aV  Ji       a^  .  2ia^       a**     ^     Ilia' 


3a    ^  2ia^       3a^  ^ 


135n 


■K- 1  a     ♦   •  •  • 


2ay      L        5a\  I     ,       2aV        .,, 


We  can  also  solve   (2U)  and  (25)   approximately  in  the   same  way  as  we   solved 
(l6)  approximately,   i.e.,  we  take  all  b's  and  c's  zero  except  b^,  b^^,  c^  and  c^. 
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The  corresponding  values  of  the   scattering  coefficient  correction  are  given  in 
Table  3.     The  expressions  for  L  ,  L,,  M    and  M-,,  necessary  for  this  calculation 
are  given  by  equations   (B9)  -  (Bi2)  inclusive, 

a  Coefficient  correction  a 


0.1 
0.2 
0.3 

0.14 

0.5 
0.6 
0.7 
0.8 
0.9 
i.C 
1.1 
1.2 
1.3 
l.U 
1.5 


0.0003010 

0.00121U 

0.002768 

0  .C0501U 

O.CO8025 

0 .01190 

0.01675 

0.0227U 

0 .03002 

0.03880 

O.CU925 

O.C6151 

O.C7571 

0 .09151 

0.1087 


Coefficient  correction 

1.6  0.1262 

1.7  0.1U27 

1.8  0.1565 

1.9  0.1656 

2.0  0.1686 

2.1  0.1650 

2.2  0.1555 

2 .3  0  .llil2 
2,k  0.1235 

2.5  0.1069 

2.6  0.08970 

2.7  0.07U31 

2.8  0.06198 

2.9  0.05U31 
3.0  0.05208 


Table  3 


The  column  headed  'Coefficient  correction •  gives  z  /2  times  the  term  which 


must  be  added  to  the  plane  wave  scattering  coefficient  to  obtain  the  scattering 
coefficient  for  a  point  source  such  that  z  »  1  and  a/ 
suraed  that  only  a  ,  a,,  b  ,  b, ,  c  and  c,  are  non-zero, 


coefficient  for  a  point  source  such  that  z  »  1  and  a/2z  «  1  when  it  is  as- 


When  a  is  small  we  can  uake  use  of  equations  CB16)  -  (B21).  The  corresponding 

values  of  b  ,  b^ ,  c  and  c,  are 
o   1   o     1 
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°  8ln^  2025n^  °  9r'^   V  >^ ) 

135n  li>n 

^  25n  -^  225»i^ 

Ttiese  results,  together  with  (18)  and  (19),  show  that  the  correction  to  the  plane 
wave  scattering  coefficient  is  approximately 

2  o  U 


2      /    8a^     ^     2aM 
o 


7,   An  upper  bound  for  the  scattered  amplitude 

In  this  section  we  give  the  general  theory  for  deriving  an  upper  bound  for  the 
scattered  amplitude  and  scattering  coefficient  in  any  scattering  problem  in  which 
the  non-radiating  problem  can  be  solved.  Since  the  non-radiating  picture  is  appro- 
priate for  low  frequencies  it  may  be  anticipated  that  the  upper  bound  will  be  most 
valuable  at  low  frequencies.  As  an  illustration  rough  upper  bounds  for  the  scattered 
■■plitude  and  scattering,  coefficient  of  a  circular  disc  in  a  normally  incident 
plane  wave  will  be  determined  and  in  the  field  due  to  a  point  source. 

The  mathematical  problem  in  general  scattering  theory  (for  sound  waves,  elec- 
tromagnetic waves  and  beams  of  atomic  particles)  may  be  presented  as: 

Find  the  solution  g  of 
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(28)  f     =     T  g 

where  f  is  known  and  T  is  a  linear  operator.' 
The  scattered  field  can  be  written  as 

[h,g]     -      [g,h] 

where  h  is  known  and  the  brackets  indicate  an  inner  product,  the  square  brackets  being 
used  to  allow  of  its  being  distinguished  from  the  inner  product  in  Section  6.  Illus- 
trations of  how  particular  problems  can  be  written  in  this  way  can  be  found  elsewhere 
(Jones  t^^ ) . 

Let  T  »  T  +  iT.  where  T  and  T.  are  real  linear  operators,  i.e.,  they  produce 

real  functions  from  real  functions.  Suppose  that  t"  ,  the  inverse  of  T  ,  exists. 

r  r 

The  finding  of  t"     corresponds  to  solving  th^  non-radiating  problems  for  the  cir- 
cular disc.     The  application  of  t"    to  (28)   supplies 

(29)  T'^h     -     (1  +  is)  g 

where  S  a  T~  T . .     Observe  that  S  is  a  real  operator* 
r     1 

Define  the  norm    ||g||  by  the  equation 

l|g||     -     [g,g*]  ^/' 
where  the  asterisk  indicates  a  complex  conjugate.     Then 

(30)  ||(l  +  iS)g||^     .      [(1  +  iS)g,   (1  -  iS)g*]  -    ||g||^     +      ||Sg||2 

provided  that  S  is  a  symmetric  operator,  i.e.,  provided  that 
[Sh^^,  h^]     .      [h^,   Sh^]    . 

7     The  analysis  is  applicable  to  vector  fields  with  suitable  interpretation  of  f ,   g 
and  T. 
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It  follows  from  (30)  that 

||g||^    <      ||(i  *  iS)g||2  . 

Hence   (29)   shows  that 

(31)  lie  11^   <    ||t;^  f||^  . 

Therefore 

(32)  I  Oi,g]  I  <  ||h||  iigji  <    ||h||  iiT;^f||. 

The  inequalities  (31)  and  (32)  provide  upper  bounds  for  the  norm  of  g  and  the 
scattered  amplitude  respectively  when  the  operator  T~  T.  is  symmetric. 

A  related  inequality  which  may  be  obtained  in  a  similar  way  is 

(33)  I  |h,g-  T^^tj  I  <  ||h||    iiT;^f||. 

A  modification  of  (32)  is  useful  in  certain  circumstances  such  as  the  calcu- 
lation of  the  scattering  coefficient.  Here  one  may  require   ^  (h,g]  (say)  and 
know  that 


(3U)  ^  [h,  T^^fl  -  0. 


The  upper  bound  given  by  (32)  is  then  usually  much  too  high.     To  overcome  this 
difficulty^ put  g  »  t"  f  +  g., where 

-  iS  T'^f     -     (1  *.±S)g^  . 
Then 

/e[h,g]  -  /^&>gj 

on  account  of  (3U)   and,   therefore. 
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OS)  i/e[h,g]  1 1 1  Ch,gj  I  <  iih||  iisT;^f  ii. 

In  general  this  is  a  better  upper  bound  for     /^[h,g]  than  that  supplied  by  (33), 
This  real  part  may  be   replaced  by  the  imaginary  part  without  affecting  the  analysis 
provided  that  the  analogue  of  (3U)  holds. 

As  an  application  of  the  above  inequalities  we  obtain  a  rough  upper  bound  for 
the  field  scattered  by  a  circular  disc  in  a  nonnaJly  incident  plane  wave.     Equations 
(1)  and  (13)  may  then  be  identified  with  (26)  and   {29)*     The  inner  product  is  de- 
fined ^1  1 

[h,g]     -     I     h(p)  g(p)  p  dp  -  J  H(y)  G(y)  y  dy 
0  o 

2  2 

where  y     ■  1  -  p     and  the  capital  letters  indicate  functions  obtained  by  the  change 

of  variable.     The  second  form  enables  us  to  make  the  identification 

1+y 


r  n    J 

1-y 


sinh  ax        , 
dx 


and 


Sp  -  -  1  1    J   sinh  a(x-y)  _  sinh  a(x^y)         ^ 


Obviously  S  is  symmetric  and  therefore  formulas  (31),   (32)  and  (35)  are  available. 
Since 

sinh  ax 


/•3A^                        smn  ax.     ^  , 

(3o)  <     a  cosh  ax. 


(37)  |t'  f I     <    -    cosh  a  sinh  oy  • 


Consequently,  (31)  implies  that 
1 

f      y|G(y)|^dy    <    -^    coEh^aR(a) 
Jo  ""  °- 
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1  2 

where      R(a)     -     a(siiih  2o  -  2a)  -  j-  (cosh  2a-l-2a  )■ 


The  scattered  amplitude  at  infinity  in  the  direction  of  the  unit  vector  n 
may  be  taken  as     e  °-  -,i  .     Therefore,  inequalities  (32)  and  (37)   show  that  the 
modulus  of  the  scatterecf  amplitude  is  never  greater  than 


(38) 


na 


cosh  a  • 


To  find  an  upper  bound  for  the  scattering  coefficient  we  observe  that  the 
scattering  coefficient  is  -2/^^,£\   and  that  ^  [l,  t"  f]     ■     0.     Thus  equation 
(3U)  is  satisfied  and  we  may  employ  the  upper  bound  given  by   (35).     Now 

fia\         A  ^  sinh  a(x+y)        sinh  a(x-y)       a   I       ,       /        x  ,       /        v 

(39)         0  <  ^  ■*•  y X  -  y         1  j-<cosh  a(x+y)  -  cosh  a(x-y)  L 


Therefore 


<  -y-  sinh  ax  sinh  oy. 


1+x 


3/  1 


sinh  ox  sinh  ay  J       - — i    ^    dt  dx  <  A-    sinh  ay, 
'l-x 


where 


1+x 


A,   -  -i-j-      sinh  ax 

5"    i  Lx 


sinh  at 


dtdx 


-  -^   ^  cosh  a  JEi(2a)  -  Ei(-2a)  >■  -e"°Ei(Ua)  +  e°Ei(-Ua)-2(Y  *  In  Ua)8inhal, 
Hence,  from  (35), 


(UC) 


*1 


^MU   ^(^h)V^. 


Seme  values  of  the  upper  bovind  for  the  scattering  coefficient,   as  given  by   (UO), 
are  shown  in  Table  U.    It  can  be  sesn  that  the  bound  gives  an  indication  of 


a 

Upper  bound 

0.1 

0.00000319 

0.2 

0.0000518 

0.3 

0.000267 

O.U 

0.000870 

0.5 

0.00220 

0.6 

O.OOU78 

0.7 

0 .0093U 

0.8 

0.0169 

0.9 

0.0291 

1.0 

0.0U79 
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a  Upper  bound 

1.1  0.0763 

1.2  0.119 

1.3  0.181 
1.U  0.271 

1.5  o.uoo 

1.6  0.586 

1.7  0.851 

1.8  3.228 

1.9  1.761 
2.0  2.515 

Table  ]^ 

The  coliunn  headed  'Upper  bound'  shows  one-half  of  the  right  hand  side  of 
(UO).  It  gives  therefore,  an  upper  bound  for  the  transmission  coefficient  of 
a  circular  aperture. 

the  order  of  magnitude  of  the  coefficient  and  is  within  10 /)  of  the  correct  value 
for  a  <  0.7  and  within  50 /fc  for  a  <  1.6.     The  estimate   (and  that  for  the  scattered 
amplitude)  could  be  improved  to  some  extent  by  replacing  the  rough  inequalities 
(36)  and   (39)  by  more  refined  ones. 

To  obtain  an  upper  bound  for  the  scattering  coefficient  when  the  incident  wave 
is  supplied  by  a  point  source  we  make  the  identification 

T'-'-  f     -     F 
where  F  is  given  by  (21),     The  scattering  coefficient    is, from  (22), 

where  the  inner  product  is  defined  in  the  same  way  as  for  the  incident  plane  wave. 
Also 
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Jt 


,    ,,     2     2,i/2 
ia(l-y  +z^) 


L    '  ^\        (i.y^z2)V2  J 


so  that  ire  can  use   (35)  with   j>^  replaced  by  ^  ,     Hence  the  scattering  coefficient 
is  never  greater  than 


(u.^,V^ 


a  ((>..^)^/^-.„) 


o 

Now,  a  use  of  (39)  gives 


,    .,     2     2,1/2 
ia(l-y  +z^)  ' 

8       e 

o        (1-y  +2^) 


SF      . 


|SF|     <     B  sinh  ay 


where 


(Ul) 


1  1+x       2         2  I2 

i°       (  sinh  ax   (        (t^l-x  )cosh  at     ^^  ^ 

3"      I  i-x    (t  *l-x^)  *UtS^  j 


o    (     .    V,  [  t  Sinh  at  ,,   ,      '2 

3n       y^  y^_^     (t  ♦l-x  )  +Ut  Zq 


Hence 


(U2) 


IISFII     <     ^    R^/2    . 


Also 


.    /,     2     2x1/2 
ia(l-y  +z   )   ' 


3       e 


3z 


o     (±-y 


,     2     2,1/2 
i-y  +zJ  ' 


*«     y 


2 

a 

2     2,2 


(1-y  >2q) 


1 

2     2J 


(1-y  +Zo) 


dy 


^  '2a^(l*2o)   *  2  ♦  \ 

^o 


?~? 
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Combining  this  result  with  (U?)  we  find  that  the  scattering  coefficient  can 
never  exceed 

'  1/2 


2a^(l+z^)  *  2  +  ^ 


BR 


1/2 


Ua^(l..^)V^  J(l.z^)V2  .  .^ 


The  actual  evaluation  of  B  from  (Ul)  is  not  easy  but  there  are  several  ways 
of  finding  an  upper  bound  for  B.     One  of  these  will  now  be  described.     We  have 


(a3). 


I6az 


3n 


sinhox 


1+x 

t  sinh  at  ...     ^     ha 

^_^       (t  ♦l-x  )  +Ut  z^  3n  z 


1+x 


sinhat 


dtdx  <  -^ —  1   sinhox   |  "'^l"'"'"  dtdx    <  A,  /z 


0    o 


l-x 


and 


1+x 


8a 


2         2 
(t  +1-X   )coshat 


1+x 


2a 


sinhax   I  'p   ' "  ^i  ^"""'g^g  dtdx  <  — ^-^      sinhax 


l-x 


coshat  dtdx 


(UU) 


where 


<  B^/z 


2a 


sinh  ax 


3n 


1+x  „ 


cosh  at  dt  dx 


-^  1(1+  -)e''^i(Ua)+(l-  -)e°Ei(-Ua)-2(cosh  a-  -^  '^)(y  +  In  Ua) 
+  2(sinh  a-  2£^)  JEi(2a)-Ei(-2a)l  -U(^~|2.  -1)  cosh  a]   . 


Thus  the  scattering  coefficient  can  never  exceed 
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The  quantity  A^(R/8a  )  '  is  tabulated  in  Table  k  and  B  {R/da^)^^^  is  given  in 
Table  5  so  that  the  upper  bound  for  the  scattering  coefficient  given  by  (U5)  can  be 
computed  for  a  given  z     without  great  difficulty. 


a 

B^(R/8a2)^/2 

a 

B^(R/8a^)^/2 

0.1 

0.0000320 

1.1 

0.0857 

0.2 

0.000261 

1.2 

0.126 

0.3 

0.000908 

1.3 

0.13U 

O.U 

0.00225 

l.U 

0.266 

0.5 

0.00it63 

1.5 

0.382 

0.6 

0,0085U 

1.6 

0.5U5 

0.7 

0.0 1U6 

1.7 

0.77U 

0.8 

0.0239 

1.3 

1.097 

0.9 

OJD375 

1.9 

1.551 

l.C 

0.0573 

Table    5 

2.0 
> 

2.191 

With  regard  to  the  actual  value  of  the  scattering  coefficient  we  may  note 
that,  if  B  were  calculated  exactly  by  (Ul)>  we  would  expect  the  upper  bound  obtained 
to  be  in  excess  of  the  correct  value  by  a  percentage  of  the   sane  order  as  that  by 
which  the  upper  bound  in  Table  U  exceeds  the  correct  value  because  the  same  method 
is  used  in  both  estimates.     The  upper  bound  given  by  (U5)  is  still  further  in  excess 
because  of  the  use  of  the  inequalities  (U3)  and   (UU).     It  is  easy  to  see  that  the 
additional  increase  is  less  than  100/(1  +  *o^^°  °^  *^®  value  given  by  (U7).     In 
particular,    if  z     >  10  the  increase  is  less  than  I /o  and  if  z     >  10  the  increase 
is  less  than  10/6.     Consequently,  it   is  suggested  that  an  estimate  of  the   scatter- 
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inp  coefficient  can  be  obtained  by  taking  the  same  percentage  of  (ii5)  that  has  to 

be  taken  of     the  right-hand  side  of  (UO)  to  obtain  the  correct  value  for  the  plane 

wave   scattering  coefficient.     This  estimate   should  not  be  in  error  by  more  than  1% 

(of  itself)  for  z     >  10  and  should  differ  from  the  correct  value  by  not  more  than 
0   —  '' 

10%   (of  itself)  for  z    >  10, 

8,       An  estimate  of  the  error  introduced  by  approximation  of  the  kernel 

Since  we  have  replaced  the  original  integral  equation  of  the  first  kind  (l)  by 
an  integral  equation  of  the  second  kind  (10 )  we  are  in  a  much  stronger  position  than 
before   for  estimating  the  error  made  by  approximate  methods.     For  example,  when  we 
solve   (10)  approximately  by  appro;:imation  to  the  kernel  as  in  Sections  5  and  6  we 
can  rewrite  the  exact  equation  so  that  one  term  involves  the  difference  between  the 
exact  and  approximate  kernels.     When  this  difference  is  not  too  large,  the  value 
of  this  term  can  be  determined  by  iteration,  and  the  exact  equation  becomes  an  inte- 
gral equation  which  can  be   solved  by  solving  a  finite  number  of  simultaneous  equations, 
The  difference  between  the  exact  solution  and  the  approximate   solution  may  therefore 
be  estimated  from  a  knowledge  of  the  magnitude  of  the  difference  between  the  exact 
and  approximate  kernels  and  a  knowledge  of  the  inverse  of  the   approximate   kernel. 
Other  processes  can  be  devised, and  the  one  that  we  employ  takes   advantage  of 
the  results  proved  in  the  preceding  section.     Attention  will  be  restricted  to  the 
problem  in  which  the  incident  wave  is  plane. 
Equation  (13)  can  be  written  as 

(U6)  (1  +  iS)G     -     F, 

Define  the  inner  product    |H,G]  by 

1 


|H,G]     =      f     H(y)G(y)dy, 
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so  that  the  scattering  coefficient  for  an  incident  plane  wave  is  -'^l^\^,fr\,    '  When 
equation  (U6)  is  solved  approximately  by  approximation  to  the  kernel  we  solve 


(U7)  (1  +  iS^)G^  -  F 

and  take  the   scattering  coefficient  to  be  -2|^[y,G  "] . 

Let  G  ■  G     +  G, .     Then  the  error  made  in  the  scattering  coefficient  by  our 
0  1 

approximation  is  -2^^,G,]  •     The  modulus  of  the  error  does  not  exceed  2  ||y  ||     jJG,  ||, 
i.e.,  2||Gj^||   /v/3. 
Now 

(1  +  iS)G^     "     (1  +  i£)   (G-Gq)  -  F  -  (1  +  iS)G^ 
or 
(U8)  (1  +  iS;G^     -     -i(S  -  S^)Gq  , 

from   (U7).     Equation  (U8)  is  of  the   same  type  as   (29)  and  ?  is  symmetric.     Hence 
an  inequality  of  the  type   (31)  holds,  namely 

IIGJI     <      ||(S-   S^)G^    II. 

Hence  the  modulus  of  the  error  made  in  the  scattering  coefficient  does  not  exceed 
2  ||(£  -  S  )G  11/  yr.  Since  all  quantities  are  known  in  this  formula  we  have  an 
explicit  expression  for  a  bound,  to  the  error, 

"I  Only  slight  modification  of  the  following  analysis  is  required  if  y  is  replaced 
by  the  general  function  H.  Thus,  effectively,  we  have  a  method  for  estimating  the 
error  in  calculating  the  scattering  coefficient  with  an  arbitrary  incident  field  and 
in  calculating  the  scattering  field  at  any  angle  from  the  solution  of  the  integral 
equation  with  the  approximate  kernel.  Note  also  that  the  method  assesses  the  accuracy 
of  the  approximate  pressure  distribution  since  it  involves  an  estimate  of  the  norm  of 
the  difference  between  the  exact  and  approximate  distributions. 
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A  reasonable  estimate  for    ||(S  -   ?-)G    ||  may  be  obtained  as  follows.     VJe  take 

„     N-i     2in+l        }  j2ni+l       .    . 

o 


c 

o 


Then,  if  H(x)     >    0  for  0  <  x  <  1, 

0  <  (S-S^)H  -  -  r    T^^^jj^       H(x)   -^5-    — ^^—    dx   . 

m-N  L  dx 

o 


/J  60 

^2m+l       .    ,  00  2s+2in+3     2s+l  2m+2 

d  sinh  ax     _    r~       a  x  ^     o „.  .      _ 


i sinh  ox  r~       a _x o 

l!m+I  X  "    ^    (i!s+i)i   (2s+2m+3)    -    2m+3 


Hence  0  <  (S  -  S  )H 
—  o 


00         2in+l  2m+2 


(U9) 


o 
1 
tl^JT    /h(x)   sinh  ox  dx. 


vhere 

N-1 


P(y)  »  sinh  oy-^    a  ^*     y     *     /  (2in+l)J 


Consider  the  special  case  in  which  N  «  2;  this  is  the   value  employed  in  Section 

5.    We  have 

■L*x  „     2  3 

2i    (  sinh  at     .+  ,  ,  ^  ,  ^° 

G       ■  -  —  X dt  +  a^x  + —       • 

o  n    J  t  0  6 

1-x 

Let  a     =  a  +ip     and  a     «  a, +ip^ ,where  a  ,   a,,p     and  p     are  real.     An  inspection  of 

the  values  obtained  reveals  that  a     and  a,    are  negative  whereas  p^  and  p^  are  posi- 
tive for  all  values  of  a  <  3.     Thus,  for  a  <  3>  £o^  <  0  ^*ien  0  <  x  <  1.     More  care 
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is  required  with  the  imaginary  part  of  C  since  it  involves  two  terras  of  opposite 
sign.  However,  it  can  be  verified  without  difficulty  that,  for  a  <  2.U,  mG^  <  0 
for  0  <  X  <  1.     Hence,  when  a  <  2,U,    (U9)   fives 


(50) 


2a 


Ol-(S-.^o)^°o     <     7?    ^2^^^^ 


and 
(51) 


2a  /3n 


0  1-(S-S,)^G,    <     ^  (g  A,  -  A3)  P(y) 


3  3 
where  now  p(y)  ■  sinh  ay  -  ay  -  a  y  /6,  and 

1 

2^3 


Ap     "         (a  X  +  a,a  x^/6)  sinh  ax  dx 


1       I  2  c 

-   ^   (a^  +  o^  +  Gj^o  /6)   cosh  a  -   (o^  +  a^  +  o^a  /2) 


2 ,^v  sinh  o 


and 


A       "  I   ^   (Po  *  ^1  *  h'^  /^^  °°^^  a  -   (pQ  ■"  Pi  *  Pl^  /^^ 


2  ,„v   sinh  a 


Therefore 


2a         .2  ^   /3n 


2    1   ^     „2" 


(S-S,)GJ|     1    ^    ^A^^(gAl-A3)'  \ 


where 


Hi  -  1  <P(y)r  dy  • 


1  /Sinh  2a 


X  (Sinn  ^a  _^).2  j(2*  |-)cosh  a-(^  +  J)sinh  a 


2         U  6 

a         a  o. 


Consequently,  when  a  <  2.U,  the  error  in  the   scattering  coefficient  made  by  our 
approximation  is  less,  in  magnitude,  than 
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(52) 


ha 


4  *{^h-  hf 


1    1 


The  values  obtained  from  this  formula  are  displayed  in  Table  6,     It  can  be  seen 
that  the  method  gives  adequate  estimates  of  the  error.     The  method  can,  of  course, 
be  adopted  to  cover  a  G^  which  changes  sign  by,   for  example,  replacing  G     by  its 
modulus  or  the  maximum  of  its  modulus  in  the  derivation  of  inequalities  similar 
to   (50)  and   (51). 


Error  bound 

0.1 

12 
0.0     5U 

0.2 

O.O^IU 

0.3 

0.0^28 

o.u 

0. 0*^38 

0.5 

0,0^2U 

0,6 

0.0000011 

0,7 

0 ,0000037 

0.8 

0.C00012 

0.9 

0.000023 

1,0 

0.000077 

1.1 

0.0002 

1.2 

0  .ooou 

Error  bound 

1.3 

0.0009 

i.l 

0.0015 

1.5 

0.0028 

1.6 

0.0050 

1.7 

0.0086 

1.0 

0.015 

1.9 

0.0  2U 

2.0 

0.037 

2.1 

0.057 

2.2 

0.C8U 

2.3 

0.12 

2.h 

0.17 

Table  6 
The  column  headed    'Error  bound'  gives  one-half  of  the  quantity  shown  in  (52). 
It  therefore  gives  a  bound  to  the  error  in  the    'Approxination '  column  of  Table  1. 
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Appendix  A 
In  this  appendix  we  consider,  firstly,  the  solution  of  the  integral  equation 

1        ^   Jo  {f>^  -  t   ) 


This  equation  may  be  rewritten  as 

7 
(A2)  f,(y)     .      f     h,(x)     '^^^  °^^--)^  dx 


j     ^(^ 


(y.x)V2 


1^^^     ■  '1^-^     -7-717^ 


(y-x)^ 


■vrtiere 


Vy 

fj^Cy)     •       [       np  f(p)dp 


and 


h^(x)     -     h(  vT  )   /2  vGT     . 


Fquation  (A2)  can  be  sol-ved  by  means  of  the  Heaviside  operator  p.     Denoting 
operational  representations  by  capital  letters,  we  obtain 

2 


or 


F^(p)     -     (n/p)^/2    e-°  /^P    H^(p) 


H^Cp)     -     (P/Ti)^/^    e°/'^P      Fi(p). 


Kence 


i/2 
/.•i\             v/\           Id       (^<./^     cosh  a(x-y)  '  , 

(A3)  h^(x)     '     -    ^     \      fi(y)  ^^_^^i/2 dy, 


80  that 


V,      ^   2     2.1/2 
cosh  a(x  -y   )   ' 


(Ah)  h(x)     -     2x      I     f(y)  y     ""^"  g^  ^"ly^ dy 
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Secondly,  we  solve  the  equation 


,  cos 


{a(xV)^] 


y  (x  -y  )  ' 


This  equation  can  be  wrttten 


y 


cos  V  a(y-x)  '    i 


f^(y,     .    J     g^(x)  ^^^^^i;;.  ^ 


where 


f^Cy)    =    f  ^  (1  -  y) 


1/2 


and 


g2(^) 


g 


This  equation  is  of  the   sane  type  as   (A2)  and  hence  its  solution  is  given  by  (A3). 
Therefore 


or 


o     A      {  °°^^  i  a(y  -X  )  ^  j 


8  U)   -   -  ^  XT     f(y)  7 


In  order  that  g(l)  »  0  it  is  necessary  that 
(A5) 


5      /                    cosh  \a{Y  -X  )  '   I 
g(x)     -    -  I  J    f(y)  y    %_Vi/ij     -^     *c  . 
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Appendix  B 

In  this  appendix  we  evaluate  the  constants  required  for  the   solution  of  the 
integral  equations   (lU)  and  (23)  by  means  of  simultaneous  equations,     ',>;e  have 

2n  /on     j2m+l         .    , 

.  a  I        2n+l     d  sinh  ax      , 

^mn     '     '(^fDU^n+Dl     j     ""  ^^7^!^    x '^• 

If  we  expand  the  integrand  in  powerp  of  x  we  see  at  once  that  A^  ■  A     .     Now 

.  a  "  j^    r  d  sinh  ax~\  /«     ,  \    /^    2n  d  "     sinh  ax   , 

^mn  ■  (W)l(^n.l)l<^    [^    -^^Jx.l'  ^'"^'M  '       ^  ~^^  ' 


/ 


(Bl) 


2a^"                    Td^"       sinh  ax" 

a         . 

(2in+l)l(2n+i)i      L^^2m           x        J^_^ 

2m+l       m-l/2,n-l/2    * 

Similarly 

A 
nm 

2a^                 r   d^"     sinh  ox"! 
(2m+i)J(2n+l}l    Ldx^"           ^       Jx-1 

a          . 
J^n+l       n-l/2,m-l/2    * 

Since  A      "A       and  A„  ,  ,„  ^  ,  .^  «  A„  i  /o  „  -i  /o> 
mn         nm  m-l/2,n-l/2         n-l/<:,m-l/^ 


o/        w  2  )     2n/^„^,x    Td  sinh  axl 

2(m-nM       -  7-F5 — \\,ir>    .^  \ ,    <    a     (2ni+l)       — 5- 

mn       (2m+l)U2n+i)I    >  ^2m  x      Jj^^^. 


2m/ ^      ,v    Pd  sinh  axl  ^^  „  j  „ 


Now 


d^     sinh  ax        ,^    ,,   ^     a^^         sinh  ax         /o„,,   ^      a^^""^       cosh  ax 


(B2)        2.^5-  °^""  °*  -   (2m) 


.  Jm         X 
dx 


,   •r—     a  Sinn  a*  /o_\f   t —       o  uumi  oa 

1  r     -T^TTT     Jm*l-2s  -   ^^-^^   ^    r^s-±)l       2m*^-28 
s«0  X  s«l  X 


I 


the  last  term  being  absent  for  m  -  0.     Hence,  when  m  ji  n, 
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mn  "   (2in+i;i(2n+l)Uin-n; 


m        2n+2e  n        2in+28 


s«0 


s-0 


(B3) 


n        2ra+2B-l 
a 


n        2n+2s-l 


W(2n.l)I  g^  2^^^^33^  .(2in.l)i  g^  1^^^-^  i  cosh  a 


the  series  starting  at  s  «  i  being  absent  vrtienever  the  sununation  would  involve  s  ■■  0, 
When  m  ■  n,   (Bi)  and  (B2)  give 

2  2m  I   m        2s  in  2s-i 

^mm-  (2m°l)U'^ml)-^  II    Wt  ^^'^  "  '  g^  rfiTETT  """^^  °/^  "  ^^    Vi/2,in-l/2- 


Also,  after  using 

d sinhox 

dx 


-   (2 


in  28-1  , 

a  cosh  ax 


nKl       2s         .   . 

a  sinh  ax 


s»l  ^  X  s-0         '       x 


we  find 


2a 


2in-l         m  2s-l 


in-1     2! 
a 


Vi/2,m-i/2"  v^ETmni  E^  r?i=T7T  °°^^  ^ "  |I  r^iir  ^^"^  *"> "  ^  Vi,m-i 


Therefore 


2]B 


mm*  (2i»+i)l 


2m  m-1   23     j^     ^ 

(2m)l(^m.l)  =^^  ^  *  L  JJ^i   gSTT  *  2^  h^"^  « 


m    2b-1 


-r^r^i^irrrs  ^*h  h°^^ 


*  2m(2m+i;  *m-l,m-l 


Repeating  the  process  we  obtain 


mm  (2m+l) I 


"m     2s 
7"   ,?  . .  J  T(2m+2)  -  Y(2s+l)V  sinh  a  -  sinh  a 


m     2s-l 

a 


■  ^   (^s-l)l  i   ^(2m+2)  -  T(2s)  \  cosh  a 


2m 


*  C2m+ljl  *oo 


rince 


-  uU  - 


(BU)   A 


oo 


d       sinh  ax    , 

X  T-    dx 

dx  X 


2   sinh  a  -  Ei(a)  +  £i(-a). 


mm 


-1 


2i\ 
( 2m+i ) I 


(B5) 


m  2s  m  2s-l 

i:      Tf?TTf  (2-2)   .  f(2s.l)l  Sinh  a  .  ^   ^f^-^^ 
S'^O  I  s«»l   ^  ' 


•  JT(2m+2)  -  Y(2s)lcosh  a  -  jJ   Ei(a)   -  Ei(-a)i 


The  evaluation  of  K  for  m  -  0  and  m  -  1  is  not  difficult  by  direct  means 

m  • 

but  the  best  method  for  general  results  appears  to  be  the  following.     We  have 

1     1+x 
„  1  (     (         sinh  at     d'^  sinh  ax    ,,     , 

'^m  "   (^m+i)i  — t —  T:2i;Tr  — x —  ^^  ^ 


^-1  l-x 


dx 


2m  -^     ^*^ 

rd 1  (     (         sinh  at     d      sinh  a(x-y)^ 

^  2m     (2m+l)J 
-dy 


3x    x-y 


dt  dx 


'-1   l-x 


i. 


y  *  0 


-dy       -'J   ■  0 


where 


2  1 

y(    )  m   si"h  a(i-y)  (  sinh  at  ,    (     sinh  a(x-y)   sinh  a(i-x) 
^^  '     1-y     j     t      ~  )     x-y         l-x        * 

-o  -1 


A  straightforward  calculation  shows  that 


-U(l-y)Y(y)   -    IeI  J2a(l-y)[  -  Ei  j-2a(l+y)|  -  Ei  J2a(y-1)     +  Ei  -120(1+7)  ^  +  cj 
.   sinh  a(l-y)-    IeI  S2fT(l-y)>  -  Ei  S-2a(l+y)|  +  Ei  i2a(y-l) 


Ei  '(J2a(l*y)  ^  ♦  2  In  ^  +   D 


cosh  a(l  y) 


where 
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C  -  Ei(Ua)  -  Ei(-Ua)  -  2Ei(2a)  +  2Ei(-2a) 
and 

D  -  Ei(Ua)  +  Ei(-Ua)  -  2y  -  2  In  ka. 
Hence 

j^  .    2  Q  r:-  (2ni)l d  ""  ^  sinh  ax   d  ^"  sinh  2  ax. 


r   " 


m   (2in+i)l  (_"  ^-,  (2r)l(^in-2r)l  ^2in-2r    5E     ^2r-l  T 


^  T—        (2in) d  "  ^"  cosh  gx  d  ^    cosh  2  ax 

^^     (2r+l)H2ni-2r-l)l  ^^2m-2r-l   3E   ^ 3E 

^  1  _  d^"  cosh  ox  p-  (2m)  I  d^°"^^~^  cosh  ox  1 


Typical  values  are 

(B6)   K^  -  J-  e"°Ei(Ua)  +  j- e*^i(-Ua)  -  (y  +  In  Ua)cosh  a, 

(B7)  \'  ^   |(a^+2a+2)e'°Ei(Ua)  +  (a^-2a+2)e°Ei(-Ua)  -  U  cosh  a(co5h  2a-l) 

-  8a  sinh  a  -  2  J  {a   +2)cosh  a  -  2a  sinh  a>  (y  +  In  ka)     , 

(B8)   Kg  -  ^  |(a^+Ua^+12a^+2Ua+2U)e"°Ei(Ua)  +  (a'*-Ua^+12a^-2Ua■^2U)e°Ei(-Ua) 

-  2  7(a^+12o  +2U)coEh  a  -  Ua(a  +6)sinh  al  (y  +  in  Ua) 

2  2  2       1 

-  8(a  +7)cosh  a  cosh  2a  -  16a(a  +7)sinh  a   +  8(5a  +7)cosh  a  . 

With  regard  to  the  evaluation  of  L  ,  no  attempt  has  been  made  to  provide  a 
general  exprepr-ion.  Instead  L  and  L,  have  been  calculated  directly.  We  find 
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(B9)       L^ 

and 

(BIO)     L, 


-|e*Ei(-Ua)  -  e"'^i(Ua)   -  2a(Y  +  In  Ua)sinh  a  +  2(sinh  2a  +  2a)sinh  a' 


"  '  SnS"  K(a-2)9^i(-^a)-a^(a+2)3"^i(Ua)-2a^(a  sinh  a  -  2  cosh  a) 


•(y  +  In  ha)  +  2(a  +2)sinh  a  sinh  2c 


+  Ua  1 


(a     •♦•  2) sinh  a  -  Ua  cosh  a 


] 


y 


Also 

(Bll)   M  -  -  [u  cosh  a  -  -  sinh  a  -  ^  sinh  2o  sinh  a  +  e"Fi(-Ua)-  e'°vi(Ua) 
0   n  I  a  c. 

>-  a 

-  2  cosh  a  (y  +  In  Ua) 
and 

(Bi2)     H.   ■  ^  2(-  +  a)sinh  a  cosh  2a  -   |i>  +  -j-jsinh  2a  sinh  a  +  2( aWnh  2a  eosha 

*—  a 

2  2     fL 

-  Lba  sinh  a  ♦  U(a  -l)cosh  a  -  (2  +  2a  -  a  )eT;i(-Ua) 

+   (a  +2a-2)e"^i(Ua)  +  J(2-a   )cosh  a  +  2a  sinh  aV  (y  +  In  Ua)      • 

It  is  helpful,  when  considering  the  behavior  for  small  a,  to  have  expansions  in 
powers  of  a  available.  The  first  terms  of  the  expansions  of  the  above  constants  are 
therefore  given  below: 

(B13)  A_t;  |a^  *  ^    a^  , 


oo 


(BlU)  A^i  ^  ^  , 

J 

(Bi5)       A^i  ;=^  fer ' 


-  U7 


(Bi6)     K  ce, 


(B17) 


.         2a 


(Bia)       L^.-^  (^i^)^ 


<Bi'>    "l  »  -  B?F  • 


(B20) 


M  ~ 


Uia 

33ir 


f  ( 


2  + 


61a 


)» 


and 
(B21) 


„  ^  Uia^  /  1  ^  53a  \ 
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Appendix  C 

In  this  appendix  we  prove  the  formula 

I       J         (isine)  J       j^(Zsin©)  sin^"^*  9  cos^^  9  d© 
o 

(Ci) 


o       ""7  ^*7 


\*  ir  (  '  J^x  <(z  -Z  cos  ©)        \ 


0 

when   /QX)  >  -  J  ,    y^(n)  >  /^(M  -  1,    /^(v)  >  -  ^^  . 

Sonine's  second  finite  integral  (Watson'-  -^)   states  that 

J   (2sin&)  I„(Zcos9)    sin"*-"-©  cos^*^©  d©  =   z^Z^     °?^"'o  )  \a.^  \  /'i 

when    ^(a)  >  -1  and    /v(P)  >  -!•     Hence  the  left-hand  side  of  (Cl)  can  be  written 

Ti/2  1  ,1  n/2 

Z  sin  ©  cos         ©  J       ,(zsin©)         J     ^(Zsin^f)!       ^(Zco9©cosgr) 

l  ^-  ?  ^o      ""  ^-  ? 

,      1 
•   sin^"         (2  cos  gr  d(2  d© 

n/2  ^^  1  n/2 

-  Z   I     sin'^"^*^  cos         (H  J     .(Zsin^l)    |     J         (zsin©)I       ,  (Zco8©cosgf) 

v+  ^  X+  J- 

.   sin  ©  cos  ©  d©dgi 

on  inverting  the  order  of  integration,     A  further  use  of  Sonine's   second  integral 
formula  gives 

z      ^Z      M       sin^-^*^g  cos^VVxCZsingr)     ^^^^\  g'^vA)/'/     ^^ 

y  ^Z   ^ Lt  cos   yi ) 

o 

aiKl  (Ci)  is  proved. 
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